Theory of superfast fronts of impact ionization in semiconductor structures 
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We present an analytical theory for impact ionization fronts in reversely biased p + -n-n + structures. 
The front propagates into a depleted n base with a velocity that exceeds the saturated drift velocity. 
The front passage generates a dense electron-hole plasma and in this way switches the structure from 
low to high conductivity. For a planar front we determine the concentration of the generated plasma, 
the maximum electric field, the front width and the voltage over the n base as functions of front 
velocity and doping of the n base. Theory takes into account that drift velocities and impact 
ionization coefficients differ between electrons and holes, and it makes quantitative predictions for 
any semiconductor material possible. 



PACS numbers: 85.30.-z, 72. 20. Ht 

I. INTRODUCTION 

Fronts of impact ionization can be excited in lay- 
ered semiconductor structures such as p + -n-n + diodes 
and p + -n-p-n + dynistorspi^i&^&i' 8 -' 9 - The front passage 
fills the structure with dense electron-hole plasma and 
hence leads to the transition of the reversely biased p-n 
junction from low-conducting to high conducting state. 
Apart from microwave TRAPATT (TRApped Plasma 
Avalanche Triggered Transit) diodes—^ triggering of im- 
pact ionization fronts has been observed in high voltage 
p + -n-n + diodes manufactured from both Si £^2j& and 
GaAsi 8 -^ In modern semiconductor electronics excita- 
tion of ionization front is a unique nonoptical method 
capable to form subnanosecond voltage ramps with kilo- 
volt amplitudes. It has found numerous pulse power 
applications i 10 : 11 

The mechanism of front propagation is based on 
avalanche multiplication of carriers by impact ionization 
and subsequent screening of the ionizing electric field due 
to the Maxwellian relaxation in the generated electron- 
hole plasma. The qualitative picture of the front passage 
is well-known (e.g., see Ref. [l2| and references therein). 
In Fig. [1] we sketch the profiles of the electric field E and 
the total carrier concentration a = n + p (here n and 
p are concentrations of electron and hole, respectively) 
in the n base of reversely biased p + -n-n + diode struc- 
ture. The front propagates into a depleted region where 
the concentration of free carriers <7o is much smaller than 
the concentration of dopants Nd- In the depleted re- 
gion the slope of electric field qNd/eeo is controlled by 
the charge of ionized donors. The ionization zone trav- 



els into the depleted region with velocity Vf that exceeds 
the saturated drift velocity v ns of electrons. This is pos- 
sible due to a small concentration of free carriers oq in 
the depleted region. The multiplication of these carriers 
starts as soon as electric field becomes sufficiently strong, 
therefore the ionization zone can propagate faster than 
the drift velocity. Generation of electron-hole pairs and 
subsequent separation of electron and hole in the elec- 
tric field form a screening region behind the ionization 
zone. Here the drift velocities remain saturated. The 
space charge in the screening region is due to excessive 
electron concentration. In the plasma layer the carrier 
concentration tx p i exceeds Nd by several orders of magni- 
tude, electric field is low and corresponds to linear Ohmic 
regime. The described propagation of ionization front 
resembles the propagation of finger-like streamers into 
pre- ionized medium (see discussion in Ref. [HI). 

The analytical theory of traveling ionization fronts in 
p + -n-n + structures has been pioneered in Ref. [2| for TRA- 
PATT diodes. This classical work is based on two crucial 
simplifications that make the theory essentially qualita- 
tive: (i) impact ionization coefficients are modeled by 
step functions and (ii) electrons and holes are assumed 
to be identical. Later on the focus mostly shifted from 
analytical studies to numerical simulations! 14 i 15 i 16 The 
demand for a quantitative analytical description remains 
strong, in particular, as ionization fronts in wide band 
materials — and also operation at electric fields above 
the band-to-band Zener breakdow n 18 ' 19 have a promis- 
ing prospective. Analytical theory is also important to 
place ionization fronts in doped semiconductor structures 
into the general context of studies on front dynamics in 
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spatially extended nonlinear systems! 20 ' 21 ! 22 

This article presents a theory of ionization front in a 
reverserly biased p + -n-n + structure. We take the asym- 
metry between electron and holes in both transport and 
impact ionization into account and assume a general form 
of impact ionization coefficients. The theory determines 
the maximum electric field in the traveling front E m , the 
voltage over the structure u, the concentration of the 
generated plasma a p i and the electric filed in plasma E p \ 
as functions of front velocity Vf and the n base doping 
N,i for self-similar propagation with constant velocity Vf. 
These results determine the instant front velocity Vf as 
a function of the applied voltage u and front position x / 
in cases when the front velocity and shape vary during 
the passage. 



where 

V±{E) s v P W±v n M (?) 

Note that for most semiconductors v~(E) < 0. 

Solving Eq. (O for p , substituting this p into the ex- 
pression d t p in Eq. (5| and integrating over x, we obtain 
the conservation of the total current density in a one- 
dimensional system 

J = q[v+(E)cj + v-{E) p] +ee d t E d x J — 0. (8) 

Here the first and the second term correspond to the 
conduction and displacement components of the current 
density J. In the following we replace Eq. ([5]) by Eq. (jHJ). 



II. THE MODEL 
A. Basic equations in drift-diffusion approximation 

We investigate how an impact ionization front passes 
through a uniformly dope n base of a reversely biased p + - 
n-n + diode structure as sketched in Fig. [I] Heavily doped 
p + and n + layers play the role of contacts and are not 
taken further into consideration. The carrier dynamics in 
the n base is described by the standard set of continuity 
equations and the Poisson equation 

d t n - d x [v n (E) ■ n] -D n d 2 x n = G(n,p, E), (1) 
d t P + d x [v p (E) -p\- D p d 2 xP = G(n,p,E), 
d x E= ^(p-n + N d ), 

where n,p and v n ,v p are electron and holes concentra- 
tions and drift velocities, respectively; E is the electric 
field strength, N d — const is concentration of donors in 
the n base, q > is the elementary charge, e and £o are 
the permittivity of the material and the absolute permit- 
tivity, respectively. We use the notation v n>p (E) > and 
take the actual direction of the carrier drift through the 
signs in equations |T]) into account. The impact ioniza- 
tion term is given by 

G(n,p,E)=a n {E)v n (E)n + a p (E)v p {E)p, (2) 

where a n>p {E) are impact ionization coefficients. 
Let us introduce new variables 

a = n + p, p = p — n. (3) 

The first variable a is the total concentration of free car- 
riers, the second one p is proportional to the space charge 
of free carriers. Neglecting diffusion that does not play 
any role on the relevant scales (e.g., see Ref. [l3l ). we 
present Eqs. (fjl as 

d t a + d x [v-(E) a + v+(E)p] =2 G, (4) 

d t p + d x [v+{E)a + v-{E)p] = 0, (5) 

8 X E = S-[p + N d ], (6) 



B. Self-similar propagation of the ionization front 

We consider the fronts moving in the same direction 
as electrons drift - so called negative fronts. We choose 
E > 0, hence electrons and negative fronts move to the 
left, cf. Fig.[T] Positive fronts in n + -p-p + structures that 
move in the same direction as hole drift, can be described 
by exchanging electrons and holes in Eqs. |T]) and replac- 
ing the donor concentration by the same acceptor 
concentration N a . For the self-similar front motion with 
velocity Vf = const we get 

<x{x,t) — <t{x + Vft), 

p{x,t) = p(x + v f t), 
E(x,t) = E(x + v f t), 

where we fixed the notation such that the fronts move 
with a positive velocity vf > in the negative x direction. 
Then in the comoving frame z = x + Vft Eqs. (HI),© an d 
© become 

d,[(v f +v-)a + v+p]=2G, (9) 
J = q [v + a + v~ p] + ee vj d z E, d z J = (10) 

d z E = ^-[p + N d ]. (11) 

eeo 

The velocity of the negative front Vf can not be smaller 
than the saturated electron velocity v ns ££L As we will see 
in the next section, a constant velocity Vf = const implies 
a time independent total current density J = const. 

C. Relation between the current density and the 

front velocity 

Using Eq. (jTTJ) to eliminate d z E from Eq. (TlU|) . we find 

d z J = 0, J = qv f N d + qj,j = v + a+[v f + v-]p. (12) 

Here j is chosen to be a particle current density whereas 
J is a charge current density. Both J and j are constants 
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in space. Typically the electric field at the left bound- 
ary Bjeft = E(x = 0) is too low for impact ionization 
[a„ ! p(£'icft) = 0] but sufficiently strong to saturate the 
drift velocities so that v n (E\ cit ) = v ns , v n (Ei cft ) = v ps 
(see Fig. [IJ. Then j is determined by small concentra- 
tions no,po <C Nd of free carriers that are present during 
the front passage in the depleted part of the structure far 
away from the ionization zone: 



j = vj a + [v f +v s \ p , 
a = n a + po, p =p ~n , v; 



(13) 



± Vp S i V ns 



Since do , po <C N d the second term in Eq. (fT2|) is negligi- 
ble. Hence J « qNdVf. This result is well known — and 
physically clear: since the front propagates with velocity 
vj into a charged medium with space charge density qNd, 
the current density qNd Vf is required to neutralize the 
space charge of the ionized donors. Hence for a planar 
front, self-similar propagation with constant velocity Vf 
corresponds to a fixed total current in the external circuit 
d t J = 0. 

Although the initial carrier density in the depleted re- 
gion is low cr <C Nd, their presence is a necessary re- 
quirement for the front to propagate with a velocity Vf 
that exceeds the electron drift velocity v ns . This fast 
propagation mode is possible because these initial carri- 
ers mulitply in an avalanche like manner as soon as the 
electric field exceeds the ionization threshold. 

In the case under study the concentrations no and po 
are not permanently nonvanishing in the medium the 
front propagates into. On the contrary, a certain mech- 
anism creates these carriers in the n base just before the 
front starts to travel, and it actually triggers this propa- 
gation. This mechanism is not universal and depends on 
the specific design and operation mode of a semiconduc- 
tor device. In microwave TRAPPAT diodes the carriers 
that remain in the structure from the previous front pas- 
sage serve as initial carriers for the next passaged In 
contrast, in high voltage diodes used as pulse sharpeners 
the time period between subsequent front passages is so 
long (typically > 100 fis) that each front passage repre- 
sents an independent event 4^ Between the pulses the 
reverse voltage is kept close to the stationary breakdown 
voltage Ub. During this waiting period that can be arbi- 
trarily long the leakage current is much smaller than q j, 
so that n base is essentially empty To trigger the front 
the applied voltage u is being rapidly increased above itf,. 
Experiments show that the front starts to travel when u 
exceeds Ub several times4^ It appears that in these de- 
vices the initial carriers are generated by field-enhanced 
ionization of deep-level electron traps i 23 ' 24 i 25 The release 
of electrons bound on these deep centers triggers the front 
and provides conditions for its superfast propagation. 2 - 
Here we focus exclusively on the stage when the ioniza- 
tion front is already traveling and refer to Refs. I24H25I for 
the detailed discussion on triggering and initial carriers 
problems. We assume that concentration of initial carri- 
ers suffices for using our density model and we treat do 



as an input parameter of our model. 



D. The final set of equations 

Eq. (fT2"l) allows to express the "space charge" p via the 
total concentration a as 



P = 



v+(E) 



J 



v f + v-(E) v f +v-(E) 



dzj = 0. (14) 



Substituting Eq. (TT4"]) in Eq. ©, we obtain the final dif- 
ferential equation for a 



(vf + v~) 2 - (v + ) 2 ^ v + 
Vf + v~ Vf + v 



■J 



2G. (15) 



The impact ionization can be expressed via variables a 
and p as 



G(a, p, E)=(3+(E) a + p-(E)p 

with impact ionization frequencies 

± = v P (E)a p (E)±v n (E) a n (E) 
P ~ 2 

Using Eq. (fT4"|) again to exclude p, we obtain 



G(a,E) = f3 cS {E,v f )a + 



■3, 



PeS(E,Vf) 



Vf + V 

f3 + [v f + v~] - f3~v^ 
Vf + v~ 



(16) 



(17) 



(18) 



Finally, we express p by a in the Poisson equation pip . 
This leads to Eq. ^ below. 

The two equations for the total carrier concentration 
a and the field E 



(vf+y-y- {v + y 

Vf+l 



T + J 

Vf + V 

2f3 cS {E,v f )a 



(19) 



2/3" 



Vf + v 



7 7-i ; q 

d 7 E = cr , 



ee Vf + v 



■J - 



(20) 



ee 



Vf 



with the known functions v ± (E), /3 ± (£') and (3 c s(E 7 Vf) 
and the constant 



j = v + a + [vf + v ] p. 



(21) 



completely describes the self-similar propagation of im- 
pact ionization front. 

Summarizing, we have first substituted the variables 
(n, p, E) by the new variables (cr, p, E) , and then we have 
expressed p by the conserved total current J(or j) and 
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a. In this way we have eliminated the third equation 
front the original set |T]) using the fact that the local 
space charge density p is no independent variable when 
the conserved total current J and the local carrier density 
a and field E are fixed. The form of the final equations 
(fT9|) and (f20|) is one of two odinary differential equations 
of first order. 



E. Special and limiting cases 

Apart from the general case v n (E) ^ v p (E), a n {E) ^ 
a p (E), it is instructive to consider the following special 
cases: 

(a) v n (E) = v p {E) = v{E), a n (E) = a p (E) = a(E), 

(b) v n (E) = v p (E) = v(E), a p (E) = 0, a n {E) = a{E), 

(c) v n (E) = v p (E) = v(E), a n (E) = 0, a p (E) = a(E). 

In case (a) electrons and holes are fully equal, in cases 
(b,c) their drift velocities are equal, but impact ioniza- 
tion takes place due to only one type of carriers. For 
example, cases (b) and (c) are simplified but reasonable 
approximations for Si and SiC,— respectively. 

Another two limiting cases provide insight in the front 
dynamics for materials where the asymmetry between 
positive and negative carriers is very strong: 

(d) v n {E) = v(E), v p (E) = 0, a n {E) = a(E), a p = 0, 

(e) v p (E) = v(E), v n (E) = 0, a p (E) = a(E), a n = 0. 

Below we refer to these special cases as to cases 
(a),(b),(c),(d) and (e), also to denote the respective 
curves in figures. 



III. GENERAL PROPERTIES OF THE 
STATIONARY FRONT PROPAGATION 



E n s,ps that are much smaller than the effective thresh- 
old of impact ionization. Consequently, in the range of 
electric fields where impact ionization sets in {a n>p =/= 0) 
the velocities of free carries do not depend on electric 
field [v+(E) = v+, v-(E) = v- in Eqs. (1191201211) ]. On 
the other hand, the generation term is negligible in the 
range of electric fields where the nonlinearity of v n , P (E) 
dependencies is essential. 

For semiconductors where the drift velocities v n {E) 
and Vp(E) depend monotonically on the field, the ap- 
proximations 



Vn,p(E) = V n 



E 



'E + E„ 



(22) 



ns,ps 



and their modifications are widely used^ The impact 
ionization coefficients a n (E) and a p (E) for electrons and 
holes are usually modeled as^ 



a n,p{E) = a„o,poexp 



E n Q,pO 

E~ 



(23) 



This approximation is known as Townsend approxima- 
tion in gas discharge physics.— The characteristic trans- 
port fields E nStPS are typically much smaller than the 
impact ionization fields E n ^^\ E ns ^ ps <C E n o tP Q. For 
example, in Si, we have E nStPS ~ 10 4 V/cm whereas 
E n o,pO ~ 10 6 V/cm.— It should be noted that none of 
these explicit approximations (f2"2")l or (f2"3")) are needed for 
our analytical results. 

Eqs. (JT9J) , ((20)) and J2l]) will be solved separately in the 
range of strong electric fields E s < E < E$ where impact 
ionization takes place and the range of moderate-to-low 
electric field E < E s where transition from high-field 
transport to low-field transport occurs and the plasma 
layer is formed. Due to the overlap between these two 
regions the obtained solutions can be sewn together pro- 
viding a consistent description. 



The main parameters of the traveling front are the 
maximum electric field E m , the width £ p of the screening 
region, the voltage u across the n base, the carrier con- 
centration CTpi and electric field E v \ in the electron-hole 
plasma behind the front (Fig. [T]). [In Fig. [2] we sketch 
the respective a(E) dependence that follows from E{x) 
and g{x) dependences shown in Fig. Q]] In this section 
we relate these parameters to the front velocity Vf, the 
doping level 2V<j and the initial concentration ctq in the 
deleted region, assuming that Vf = const and front prop- 
agation is self-similar. Note that according to Eq. (fT2|) . 
the velocity Vf can be expressed by the current density 
J = qN d v f . 



A. Ordering of scales 

In semiconductors the drift velocities v n>p (E) typically 
saturate in electric fields above the characteristic fields 



B. Equations in the high field region 



For E > E ns ^ ps we assume that carriers drift with con- 
stant velocities v ntP {E) = v nStPS . Then Eqs. (|19I20I21|) 
become 



d z a = \p cS (E : Vf)a, 

d z E = b — ccr, 

j = v+a + (v f +v~)p, 



where 



A = 



2(v f + v-) 



{v f +vs) 2 -(vf) 



ee 



Vf + v s J ee 



qN d 



q vj ± v ps ± v ns 

c = , vf - 



(24) 
(25) 
(26) 



(27) 



(28) 



eso Vf + v s 
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The coefficients A, b and c are constants. When deriving 
Eq. (f2"4")l from Eq. (JTSJ) we neglect the second term on 
the right-hand side. This is justified because the order 
of magnitude value of this term is f}~ (vf /vf)ao which 
is much smaller than the first term. We also take into 
account that j/(vf + v~) *~ (v s /vf)ao <C Nd in the ex- 
pression for b. 



C. Carrier concentration just behind the ionization 

zone 

At the point Xf where electric field reaches the maxi- 
mum value -E m the slope of electric field profile is equal to 
zero (see Fig. [T]). In the comoving frame xt corresponds 
to the point z m = Xf + Vf t that does not vary in time. 
It follows from d z E = and Eq. ((25]) that 



er m = a(z m ) 



(29) 



V f + V S N , J „ Vf + Vs AT 

T JV d + —r « t l\d- 



Dividing Eq. l|2"4")) by Eq. ([23)1 we exclude z, and then 
employ (125]), This yields 



dcr = — (3 c s(E : vj) dE. 



(30) 



Let us denote as cr* the concentration cr which is reached 
just behind the ionization zone (Fig. [1]). Then the in- 
tegrals of the left-hand side of Eq. (f30|) from ao to cr m 
and from cr* to cr m are both equal to the integral of the 
right-hand side from to E m (see Fig. [2]): 



-da = 



-da = 



(31) 



p eS (E,v f )dE. 



Taking integrals over a, we find the relation between the 
concentration behind the front a*, front velocity Vf/v s , 
the doping level Nd and initial concentration <7o' 



* t a 1 a 

a = c m m h a ~ a m in — . 

Co Co 



(32) 



The concentration a* does not depend on f3 D g and thus on 
the specific form of a^ p (E). Equation (|29|) generalizes 
equation (29) in Ref. 0. 

In Fig. [3] we show a* and a m as functions of Vf /vf for 
different values of ao/Nd- Solid lines 1, 2, 3 correspond 
to the symmetric case vj = 0. These dependencies are 
valid for all three special cases (a,b,c) (see Sec. 2E) if 
vf is replaced by v s . We see that a* is approximately 
10 times larger than tr m (line 4) whereas a m exceeds ao 
by several orders of magnitude. Dashed lines Id, 2d, 
3d and dotted lines le, 2e, 3e show a* for two extreme 
asymmetric cases v^ /vf = =pl that correspond to im- 
mobile holes [case (d), v ps — 0] and immobile electrons 



[case (e), v ns — 0], respectively. In the case (d) of im- 
mobile holes a* — > when Vf/vf — > 1 (although the 



part Vf < 2v]. 



is unphysical). In contrast, in the 



case (e) of immobile electrons a* practically does not 
depend on »j for vj/vj < 1. This plateau exists for 
v~ > and corresponds to the interval of front velocities 
v ns < Vf < v ps (recall that for negative front vt > v ns ). 
Even for the limiting cases vj /vf = =pl the effect of 
transport asymmetry is small for Vf/vj > 10. Actual 
values of [v^ /vf\ are much smaller: at room tempera- 
tures we have vj /vf ~ —0.1,-0.05,-0.3 for Si, GaAs 
and SiC, respectively^ 

According to Eq. (f2"o| the space charge of free carriers 
behind the ionization zone p* is proportional to cr*: 



P = 



J 



v f 



"I 



"I 



(33) 



Electron and holes concentrations n*,p* = (cr* =p p*)/2 
are recovered as 



Vf + v ps 



2 Vf + v ps - v ns 

Vf ~ Vns 
2 Vf + Vps - V ns 



(34) 



The appearance of negative space charge p* — p* — n* < 
is a combined effect of spatially inhomogeneous ionization 
and separation of electrons and holes in strong electric 
field. 



D. Maximum electric field 

The maximum field E is determined by integrating Eq. 
(E01): 



A \ (Tq C m 



(35) 



Substituting the expressions for (3 c g(E), A, b and cr m , we 
obtain the explicit formula 



v f 



a n (E) dE- 



(36) 



Vf + v s + vj J 



qN d 



In 



a p (E) dE = 

Vf + v~ N d 
vf c 



Straightforward numerical integration of Eq. (136[) for 
given dependencies a n , P (E) makes it possible to deter- 
mine E m as a function of Vf and ao /Nd for any semicon- 
ductor material. 

For the special cases (a,b,c,d,e) (see Sec. 2E) and the 
Townsend's dependence a(E) = ciq exp(— Eq/E) Eq. 



6 



([36]) yields 



L 



exp \ — )dy= , (37) 



E* = 



E* = 



E* = 



7 




2 2 
«/-«« 


Of) 




b 






v s 


b 


Vf T v s 


ao 


Vs 



ln[^ 

V s <T 



So 

- 1 



for case (a), (38) 



l n [^]-l 

V s CT 



In 



for cases (b, c), (39) 
2v f TVs N d \ _ 
v s a J 

for cases (d, e). (40) 



The dimensional coefficient b/ao in E* has a meaning of 
change of the electric field on the length of impact ion- 
ization Qfg -1 for the slope d z E = b which is determined by 
the doping level Nd- For realistic parameters b/ao *C Eq. 
Comparing Eq. (|39[) with Eq. ([40]) we see that asymmetry 
of transport properties has logarithmically weak effect on 
E m . 

We solve Eq. (|3T|) numerically and show E m /E* as 
a function of E /E* in Fig. [4] Together with Eqs. 
(|38I39I40[) this dependence makes it possible to deter- 
mine £" m for given values of Vf/v s , ao/Nd, ao and Eo 
for the special cases (a-e). [Note that E m /E* — > 1 
when Eq/E* — > 0. This limit corresponds to ultrastrong 
electric field when a(E) — > ao-] E m increases with Eq 
because the effective threshold of impact ionization be- 
comes higher. 

In Fig. [5]we show E m as a function of Vf/v s for dif- 
ferent values of ao/Nd- Thick solid lines la, 2a, 3a cor- 
respond to the symmetric case (a) [Eq. (|3"5|) ]. thin solid 
lines lb, 2b, 3b and dashed lines lc, 2c, 3c correspond 
to the cases of "monopolar ionization" (b) and (c) [Eq. 
(f3"9"| ] , respectively. E m is the smallest for the symmetric 
case a n = a p because both types of carriers are involved 
in impact ionization. Curves lc, 2c, 3c are "discontinu- 
ous": E m tends to a finite value when Vf /v s — > 1 whereas 
we expect E m — for Vf/v s = 1 as it is for the curves la, 
2a, 3a and lb, 2b, 3b. This feature results from the as- 
sumption a n — and disappears for arbitrary small but 
nonzero value of a n . Physically it means that impact 
ionization by electrons that move parallel to the front 
is important at low front velocities Vf ~ v s even in the 
case a n <C a p . The curves for two types of "monopo- 
lar ionization" (b) and (c) become close with increase of 
Vf. The dependencies E(vf) calculated for fully immo- 
bile holes or electrons [cases (d) and (e), Eq. (140]) . the 
respective curves not shown] turn out to be very close 
to the cases of monopolar ionization (b) and (c) [Eq. 
(f3"9")) ], curves lb, 2b, 3b and lc, 2c, 3c, respectively]: 
the difference is smaller than 2% in the whole interval 
1 < Vf/v s < 100. It means that E m is much more influ- 
enced by the asymmetry of impact ionization coefficients 
than by the asymmetry of drift velocities. We also see 



that _E m and hence the effective width of ionization zone 
If decrease with ao- This observation explains why the 
concentration a* (see Fig. [3]) also decreases with ct - In 
Fig. [5] we choose b/(a E ) = 0.0002. For Si this value 
corresponds to N d — 10 14 cm~ 3 which is a typical dop- 
ing level for high voltage Si structures'' 4 -'^ Since in Si 
impact ionization by electrons dominates and v ns ~ v ps 
f Ref . I27T] . curves lb, 2b, 3b provide a good approximation 
for this material. 

Fig. [5] shows the dependence E m (vf) for different val- 
ues of b/(aoEo) and the fixed value of ao/Nd- E m in- 
creases with Nd due to the decrease of the effective width 
of ionization zone If and decreases with ao due to the 
more efficient impact ionization. 

The E m (vf) dependencies obtained for the case of sym- 
metric ionization a n — a p and ionization by electrons 
(curves la, 2a, 3a and lb, 2b, 3b in Figs. I5I6|) can be 
fitted by the squareroot function 



Em(vf) - E th ~ E y (vf/v s ) 



1, 



(41) 



where Eth « 0.15. ..0.2 Eq plays the role of the effective 
threshold of impact ionization. This fit is quantitatively 
accurate for v / /v s > 2 and remains qualitatively correct 
for 1 < Vf/v s < 2. Straightforward examination of Eq. 

shows that the squareroot dependence corresponds 
to the piece- wise linear approximation of the impact ion- 
ization coefficient a(E) 



a{E) ~ (E - E th ) Q(E - E th ), 



(42) 



where Q(E) is the step function. The function 
(|4"2")l approximates the Townsend's dependence a{E) = 
Eo exp(—Eo/E) reasonably well in the most important 
interval of electric fields 0.3...0.7 E (see inset to Fig. 
[5]). [Previously the approximation (|4"2"|) has been dis- 
cussed in the theory of finger-like streamers in Ref. [30L ] 
It is remarkable that the dynamics of ionization fronts 
reveals the existence of effective threshold electric field 
E t h ~ 0.2 Eq in spite of the absence of any kind of cut- 
off at low electric fields in the Townsend's dependence 
itself. In particular, the squareroot dependence (|41|) im- 
plies that Vf ~ if 2 if we define If * (E m — E t h)/b. 



E. The width of the screening region 

The screening region is situated just behind the ion- 
ization zone (Fig. [T|). Here the electric field is insufficient 
for impact ionization, but the drift velocities remain sat- 
urated. According to Eqs. (|19l2ip in this interval of 
electric fields the concentration a and space charge p are 
conserved and keep values a* and p* determined by Eqs. 
(|32I33[) . The slope of electric field in this region is deter- 
mined as 



d,E 



q 

ee 



(p* + N d ). 
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Taking into account Eqs. (|29I32[) we find the ratio be- 
tween the slope \d z E\ in the screening region and the 
slope b in the depleted region 



\d z E\ 



Vf +V S 



tJ 

N~d 



ln^ 

ctq 



(43) 



The width i p of the screening region can be evaluated as 



E n 



\d z E\ b[\n(a*/a ) - 1] ' 



(44) 



Calculating l p in this way we include in it a part of ioniza- 
tion zone (see Fig. [T]) and hence somewhat overestimate 
the actual width of the screening region. On the other 
hand, Eq. (|4"4")) gives an idea of the effective front width 
since it accounts for the two regions of most dramatic 
change in concentration and electric field: the region of 
most rapid increase of concentration from a m to the fi- 
nal value a* (approximately by one order of magnitude) 
and the region of steep drop of electric field from E m to 
E p i -C E m . 

We show \d z E\/b as a function ofvf/v£ in Fig. [7] Solid 
lines correspond to the symmetric case vj = for differ- 
ent values of ao- These curves correspond also to the spe- 
cial cases (a,b,c) if vj~ is replaced by v s . The electric field 
profile in the screening region is approximately 10. ..20 
times steeper than in depleted region. The slope \d z E\/b 
has weak logarithmic dependence on dj and decreases 
with ao/Nd- Dashed lines and dotted lines represent 
\d z E\/b for two limiting asymmetric cases v~ /vj = —1 
[case (d), immobile holes] and v~/v+ = 1 [case (e), im- 
mobile electrons]. Similar to a* , \d z E\/b has a weak de- 
pendence on the asymmetry of saturated drift velocities 
(v~ =/= 0) for sufficiently large Vf. 



F. Transition from high- field to low-field region 

In the region where electric field is insufficient for im- 
pact ionization it follows from Eq. (TT5|) that 



[v f +v~(E)] 2 ~[v+(E)] 2 v+(E) 



v f +v-(E) 



j = const, (45) 



Vf + v 

where the second small term is negligible. Employing 
conservation of this quantity, we find explicit expressions 
for a and p: 



a(E) 



P{E) 



[v 2 +v-(E)] 2 -[v+(E)} 2 vf+v s 
v+(E) [v 2 f+v-] 2 -[v + (E)} 2 



■a*, (46) 



Vf+v7 [v 2 +v-{E)f-[v+{E)f 

For the special cases (a,b,c) this yields 

v 2 ' 
a(E) = - T l 

v f 



v 2 -v 2 (E) 



(47) 



P(E) 



v(E) v 2 -v 2 



Vf Vf - v 2 {E) 



If v{E) is monotonic, then according to Eq. (|47|) \p{E)\ 
monotonically decreases with decrease of E. The transi- 
tion to neutral plasma occurs when \p(E) | reaches Nd and 
E reaches a certain constant asymptotic value. In semi- 
conductors with nonmonotonic v n {E) dependence (e.g., 
GaAs) that has maximum at E < E s , we expect that 
\p(E)\ > \p*\ near E, but the transition to plasma at 
lower electric fields occurs in the same way as for mono- 
tonic v(E). 



G. Parameters of the plasma region 

Plasma concentration and electric field in plasma are 
denoted as a p \ and E p \, respectively (Fig. [J). Generally, 
these parameters are determined by Eqs. (I21I45|) together 
with the neutrality condition p p \ = —Nd- 

For the general asymmetric case v~(E) ^ we ap- 
proximate the drift velocities in plasma by the Ohm low 
v n (E) — p n E, Vp(E) — p p E (note that for the approxi- 
mations (|22J| p n<p = v nSt p S /E nS:PS ). This yields 



Vf + u_ s 



(v f +p-E pl ) 2 -(p+E pl f 
Vf + p-Epi 



j = p + a pl E p i + (v f + p E p i)(-N d ), 



(48) 



7 pi . 



(49) 



± Pp i p r 

P" o 



Expressing E p \ via a p \ in Eq. (149p , neglecting j and sub- 
stituting in Eq. Ij48|) we find explicit formulas for a p i 



a pi = -A a* 




4iV d / N d 



p_ 



A 



E, 



pi 



E+ 



V f/ V s 



Aa* \Aa* 

Vf (v f +v7) 

v+ 



(50) 



(fpi/Nd-p I p + 
Expansion over Nd/Aa* leads to 



El 



dpi = A er* 



Aa* 



» Aa* 



(51) 



(52) 



Vf (v f + v7) 



Electron and hole concentrations in plasma region are 
recovered as 



N d 



P P i 



N d 



(53) 



It can be easily derived from the continuity of electron 
and hole flows that n* > n p \ and p* < n p \ [see Eq. (J34|)]. 
In contrast, the relation between a* and a p \ is not uni- 
versal: A < 1 when v~ < 0, but A > 1 for v~ > and 
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Vf > [{v+) 2 — (v s ) 2 ] /(v s ) 2 . Therefore according to Eq. 

for any Vf. However, for 



(jpi < a* when 



sufficiently large Vf we find that a p \ > a* when v ns < v ps . 
Next, A -> for Vf — > v ns and hence according to Eqs. 
(|5QI53p (Tpi — > iVd, u p i — > A^, _p p i — > 0. This is consistent 
with condition t>/ > tVts- 

For the special cases (a,b,c) when electron and hole 
drift velocities are equal it is convenient to use the ap- 
proximation v(E) — v s E/(E + E s ) in the whole range of 
electric fields including plasma. After similar derivations 
we find 



account that according to Eq. (|5Tj) in symmetric case 
Ep\ ~ Vf/a p i : E p i decreases and increases when the 
increase of a p \ with Vf is super linear and sublincar, re- 
spectively. Hence "plateau" on the dependence <J p i(vf) 
corresponds to minimum on E m (vf). 



H. Voltage over the structure 



The voltage over the structure is given by the integral 



Cpl 



2vj 



1 + \ 



2v) 



N, 



vi a* 



E p i = E s 



D 



D 



2 2 

2v1 



\ 



2v) 



(54) 



(55) 



In 21. 
ct 



Expanding over N^/a*, we obtain simplified dependen- 
cies that are valid when a* 3> Nd 



J f 



J f 



Epi — E s 



vj a* 



Inf--1 



<r*,(56) 



J f 



CO 



The concentration cr p i tends to a* with increase of 
and, similar to a*, shows quasilinear dependence on vj 
for Vf v s . In contrast, the dependence of E p i on v/ is 
very weak. 

The concentration cr p i and the electric field E p i are 
shown in Figs. [5]and[ni respectively. In panels (a) these 
dependencies are shown for different values of ao/Nd- 
Solid curves 1, 2, 3 correspond to the symmetric case 
v~ = 0. Dashed curves Id, 2d, 3d correspond to the 
limiting case (d) of immobile holes v~ /v+ = —1. Dot- 
ted curves le, 2e, 3e correspond to the opposite limiting 
case (e) of immobile electrons v~ /v+ — 1. For fast fronts 
cfpi increases with Vf linearly, whereas the electric field 
is close to Epi w 0.1 E+ and weakly decreases with Vf. 
In panels (b) er p i and E p i are shown for different val- 
ues of v~ /vj and We see that the asymmetry 
in high-field transport is much more important than the 
asymmetry in the low-field transport. For v~ > (curves 
4,5,6,7 in Fig.[5](b)) the dependence a p i(vf) has a kind of 
plateau and the corresponding dependence E p \(vf) has 
minimum. This occurs in the interval of front velocities 
Vns < Vf < v ps . The plateau on a p i(vf) is caused by 
the peculiarity of a*(vf) dependence that has been dis- 
cussed above in Sec. 3C (see also Fig. [3]). The nonmono- 
tonic behaviour of E v i(vf) becomes clear if we take into 



E(x)da 



(57) 



We approximate the actual profile E{z) by a piece-wise 
linear profile A-B-C-D shown in Fig. [TU] neglecting the 
voltage drop over the plasma region where the electric 
field is low. For such profile <j(z) — cto and d z E = b — cao 
on part A-B, E = E m and a increases from oq to 
cr* in arbitrary way on part B-C, and o(z) = a* and 
d z E — b — co* on part C-D. Integration over this pro- 
file gives an upper bound for the actual voltage. On the 
(cr, E) plane this profile corresponds to the rectangular 
a(E) dependence (dashed line A-B-C-D in Fig.©. The 
integral over z can be replaced either by the integral over 
electric field E or over concentration a since according to 
Eqs. P4T23)) dz = dE/(b-ca) = da/\(3 cS a. Employ- 
ing integration over E for branches A-B and C-D and 
integration over a for branch B-C , we approximate the 
integral ([57]) as 



EdE 



This yields 
1 

U ~ 2c" 



cer 



- ^lcft 



Err, da 



EdE 



A /3cff(#m) o J Em b- ca* ' 
(58) 



El 



(59) 



In- 



A/3 c ff(#m) 00 

Here the first term corresponds to the contribution of 
the inclined parts A-B and C-D of the field profile. 
The second term corresponds to the horizontal part B- 
C and hence approximates the contribution of the non- 
linear part of the profile near its maximum. The ratio 
of the first term to the second one can be estimated 
as v s a(E m ) t p , where v s a(E m ) is the frequency of im- 
pact ionization, r p = £ p /vf is the time the front takes 
to move over its own width l p and it is assumed that 
v n = v p , a n = ct p , Vf 3> v s - [For this estimate we 
assume E m 3> £q c ft and employ Eq. ([44|l .] Typically 
v s a(E m ) t p 3> 1 and hence the first term in Eq. (|59"|) 
dominates over the second one. It means that the non- 
linear part of the profile E(z) which is located near its 
maximum is small, and the typical shape of E(z) is close 
to triangle. 
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Taking into account that £ieft — E m — bxf, we can 
represent u as 



bx) 

u ~ E m xj — 



E„ 



2 2&[ln(cr*/ cr o) - 1] 



(60) 



In — . 



Since the dependencies cr*(vf) and E(vf) are already de- 
termined [ see Eqs. (|32I36|) ]. Eq. (|60p gives u as a function 
of front position x f and velocity v / . 



IV. ULTRAFAST FRONTS (v f > v+) 



Then it follows from Eqs. ([6l| |62 | that 



a ee {E m - E th ) a a Q £E Q (E m - E th ) 

cr m In — w . 

q cr m q 

(66) 

Predictions for Vf given by Eqs. (I64[) and (|65|) differ 
in the definition of r. However, the relative difference 
(r — t)/t = ln(cr*/a- m )/ ln(cr*/cro) does not exceed fO % 
(see Fig. [3]). Next, according to Eq. ct p i is propor- 
tional to E m whereas more accurate Eq. (|66]) predicts 
proportionality to the difference between E m and E t h- 
Therefore Eq. (|6"4"j) overestimates er p i. Still it gives cor- 
rect order of a v \ since in practice Eth., E m and \E t h — E m \ 
are of the same order of magnitude. 



The front velocity vj is often much higher than 
W 8< ?i 4 i?i9i7i9i? In the respective limiting case Vf/vf 3> 1 
the effect of transport asymmetry vanishes. As it follows 
from Eqs. (|29I32I52[) the concentrations a m and <r* can 
be presented as 

crm = -T- Nd, a* = ^-N d \n—, (61) 
vj vj (Jo 

The plasma concentration and electric field in plasma are 
given by 



a* » a*, E p i 



El 



ln(tr*/<7 ) 



(62) 



The maximum field is determined by [see Eq. ([33)1 ] 



[l>„ s a„(.E) -I- v ps a p (£)] dE = 



(63) 



gjVd«f ln v f N d 



ee 



vf a 



Let us compare these predictions with the results of an 
elementary model suggested for planar ionization front 
in diode structures in Ref. [13 on the basis of the ideas 
developed for finger-like streamers in Ref. I3flf32l . In Ref. 
IT7I it is assumed v n — v p , Vf ^> v s , a p — and a n (E) = 
ccq Q(E — E t h), where Q(E) is the step function. Under 
these assumptions the order of magnitude values of Vf 
and dpi has been evaluated as 

vj = — , t = In — , if = , (64) 



t v s a a <to " qN d /ee 

ao£EoE m 



where r is the time it takes for the front to pass over the 
width of ionization zone if. 

For a n (E) = ao Q(E ~ E t h), ot p {E) = we obtain from 
Eq. © 



v f 



4r 

T ' 



1 , (Tm 

In — 

v s ao (To 



(65) 



V. NONSTATIONARY PROPAGATION 

A. The adiabatic condition 

The relations between vf, E m , a p \ and E p \ obtained 
for J = const still hold when J varies in time providing 
that this variation is slow in comparison with inner re- 
laxation times of the traveling front. These times are the 
Maxwellian relaxation time in plasma behind the front 
tm and the time r p = £ p /iif it the front takes to move 
over the width of the screening region i p . Indeed, any 
change of the electric field (and hence the current density) 
in the n base originates from changes of electric charges 
in the highly doped p + and n + layers that serve as effec- 
tive electrodes (see Fig. [T]). Further transfer of electric 
charge into the n base occurs through the plasma layer 
and is controlled by tm • Redistribution of charges in the 
traveling screening region at the plasma edge takes time 
T p . Thus the times tm an d t p characterize how fast the 
front relaxes to the steady profile that corresponds to the 
instant value of the current density. Employing Eq. (|44p . 
assuming for simplicity v s — v ns — v ps , /i = fi n = [i p and 
taking into account v s — /J-E s , we obtain 



IE. 

v f 



TM 



Em CTpl 



qH<T p \ 



(67) 



We see that r p is much larger than tm and hence it is 
the time r p that eventually controls the relaxation of the 
front profile. Consequently, the adiabatic condition for 
the variation of current J can be presented as 



d(ln J) 
dt 



< I 



(68) 



Below we show that this condition is typically met for 
the realistic operation mode of high voltage diodes used 
as switches in pulse power applications. 



B. Coupling to the external circuit 

In practice the device is connected to the voltage source 
U(t) via a load resistance R. The current density J and 
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the voltage over the structure u are related via Kirch- 
hoff's equation 



u(t) + RSJ{t) = U(t) . 



(69) 



In high voltage diodes used in pulse power applications 
the front passage switches the structure from the noncon- 
ducting state to the conducting state. At the moment t = 
t when the front starts to travel u ~ U(t Q ) and J « 0. 
The switching time is determined as At = W/ (vf), where 
W is the n base width and (vf) is the mean value of front 
velocity (generally, Vf increases during the front passage). 
The device resistivity after switching is negligible in com- 
parison with the load resistance R. Hence u(to + At) w 
and J (t + At) = U(t + At)/(RS) « U{t + At)/(RS), 
where we take into account that variation of U within the 
time period At is small. Then we estimate the relative 
variation of the current density as 



d(lnJ) J (to + At) - J (to) 



dt 



AtJ(t + A t ) 



1 

At 



W 



Substituting this estimate to Eq. (f6"5| . we present the 
adiabatic condition as 



( v f) 1 p ^ <-p ^ 1 



W 



If 

w 



(70) 



Eq. ([70)) states that the inner dynamics of the traveling 
front and the outer dynamics that is controlled by the 
external circuit can be separated if the effective front 
width l p is much smaller than the size of the system W . 
Using Eq. flU} we present ^TDJ) as 



1 



bW ln(a*/a ) - 1 < L 

According to Eq. (j43]) and Fig. fTJ) the second term in (|7T|) 
has numerical value in the range 0.05. ..0.1. Therefore it is 
necessary that E m /(bW) ~ 1. The adiabatic conditions 
(I70I71P are met or nearly met for high- voltage sharpening 
diodes where W ~ 100.. . 300 ^m and £ p ~ 10.. .20 fim but 
are not likely to be met for much smaller TRAPATT 
diodes. 

In conclusion, the relations between Vf, E m , a p \ ob- 
tained for the self-similar propagation mode can be used 
in general case to relate the instant values of these pa- 
rameters for sufficiently large structures and fast fronts. 
In this case the voltage u(vf, Xf) given by Eq. (|60|) can be 
substituted to the Kirchoff's equation (f6"9"| . Then equa- 
tion dxf/dt = — Vf (recall the Vf > for the front trav- 
eling in the negative x direction) together with Eq. (|S"9"|) 
represent a set of ordinary differential equations that de- 
scribe the front propagation with account taken for the 
external circuit. 



VI. SUMMARY 

Basic parameters of plane impact ionization fronts in 
reversely biased p + -n-n + structure (Fig. [1} are deter- 
mined by current density J and concentration of initial 



carriers ao (regime parameters) , doping of the n base Nd 
(structure parameter) and such material parameters as 
saturated drift velocities v ns and v ps , low field mobilities 
H n and fj, p and electron and hole impact ionization coef- 
ficients a n (E) and a p (E). The front velocity is given by 
Vf « J/qNd [Eq. I|12p]. The concentration of generated 
plasma a p \ and electric field in plasma E p \ determined by 
Eqs. (|29l32l50p do not depend on impact ionization coeffi- 
cients a n ^ p (E). Concentration a v \ weakly decreases with 
initial carrier concentration [Fig. [S^a)]. For moderate 
front velocities Vf < 5(v ns + v ps ) the concentration a p i 
and field E p \ are sensitive to the ratio v ns /v ps , whereas 
the asymmetry in low-field transport /i ra //i p ^ 1 has very 
little effect [Fig. [3(b) and Fig. [3(b)]. For higher front ve- 
locities CTpi and E p i do not depend on v ns /v ps and ^ n /n P 
[Eq. (|62[) ]: a increases with Vf quaislinearly whereas E p i 
weakly decreases. 

General dependence of maximum electric field £" m on 
Vf is given by Eq. (|36p . Due to strong nonlinearity of im- 
pact ionization coefficients a n p (E) often only one type 
of carriers contributes to ionization. In this case the de- 
pendence E m (vf) can be determined in a simple form for 
the Townsend's approximation a(E) = ao exp(—Eo/E) 
and symmetric transport v ns — v ps [see Eq. (|3~T|) and 
Fig. 2] . We reveal the existence of the effective thresh- 
old of impact ionization E t h ~ 0.2 Eq (Figs. [5] and [6]) 
and the squareroot character of the E m (vf) dependence 
[Eq. (gTJ]. Eq. (gTJ) implies that v f ~ if 2 , where If is 
the effective width of ionization zone. The squereroot 
dependence fails for slow fronts when a n <C a p . E nl 
increases with Eq and Nd (Fig. [S]) and decreases with 
ao and Qo (Figs. I4I5I6P - The width i p of the screening 
region where electric field falls from E m to E p \ weakly 
depends on vf and ao [Eq. 03] • The slope of the electric 
field in the screening region is about 10... 20 times larger 
than the slope qNd/seo in the depleted n base the front 
propagates to (Fig. [7]). 

The voltage over the structure u is determined by the 
front velocity Vf and the front position Xf fEq. I59p . The 
profile of electric field E(z) is essentially triangular since 
the nonlinear part near its maximum E = E m is small. 

In the case when the current density J varies in time, 
the front velocity Vf and the front profile E(z) are non- 
stationary. The largest inner relaxation time r p = i p /vf 
is the time it takes for the front to travel over the width 
of the screening region l p [Eq. {67}]. The relations be- 
tween basic front parameters E, £ p , a v \ and E p \ obtained 
for J = const remain valid if temporal variation of J is 
slow with respect to r p [Eq. (f6"5]) ]. For the actual case of 
the device connected in series with an external load this 
adiabatic condition can be presented as i p /W <C 1 
[Eq. I|70p]: inner and outer dynamics can be separated if 
ionization front is thin with respect to the n base width 
W . This condition is met or nearly met for high- voltage 
sharpening diodes. 

For very strong electric fields E > Eq the direct band- 
to-band tunneling (Zener breakdown) must be taken into 
account. Numerical simulations show that in presence of 
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this ionization mechanism the character of front prop- 
agation substantially changes*^ The respective fronts 
have been called tunneling-assisted impact ionization 
fronts J£ Recently the dynamic avalanche breakdown of 
high voltage diodes with stationary breakdown voltage 
Ub ~ 1.5 kV at extremely high voltage about 10 kV 
that corresponds to electric fields above the threshold 
of Zener breakdown has been observed experimentally. 19 
The analytical theory of tunneling-assisted impact ion- 
ization fronts will be reported separately. 
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FIG. 1: Sketch of the electric field E and total concentration of free carriers concentrations a = n + p (lower panel) in the p + -n- 
p + structure during the passage of the ionization front. The field E s corresponds to the transition from linear low- field transport 
to saturated drift velocities. Coordinates x and z — x + Vf t correspond to stationary and comoving frames, respectively. Note 
relations cro <C Na, a* , cr p i and between the initial concentration <ro in the depleted region, doping Nd and plasma concentration 
CTpi. The relation a* > a p \ generally holds only for v n3 > v ps and can be broken for v„ s < v ps . 




FIG. 2: Dependence of total carrier concentration on electric field u(E) in the travelling ionization front. See notations and 
comments to Fig. [1] Path A-B-C-D corresponds to piece-wise linear approximation of the field profile shown in Fig. 1101 
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FIG. 3: Concentration a* just behind the ionization zone as a function of Vf/v^ for different values of ao/Nd- Thick solid lines 
1,2,3 correspond to symmetric case vj" = 0. Dotted lines Id, 2d, 3d and dashed lines le,2e,3e correspond to the two limiting 
asymmetric cases /vt = —1 [immobile holes, case (d)] and vj jvt = 1 [immobile electrons, case(e)], respectively. Curves of 
1st, 2nd and 3rd series correspond to ao/Nd — 10~ 3 , 10 -4 , 10~ 5 , respectively. Thin solid line 4 shows concentration er m at the 
point of maximum electric field for the symmetric case vj = 0. 




FIG. 4: Maximum electric field E m as a function of Eo according to Eq. (|37[) . Both £ m and Eo are normalized by 
E*(v f /v s ,ao/N d ,a ). Note that E m /E* -» 1 when E /E* 0. 
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FIG. 5: Maximum electric field _E m as a function of Vf/v B for different values of ao/Nd according to Eq. (|3T[) . Thick solid curves 
la,2a,3a correspond to the symmetric case (a) a n (E) = a p (E), v ns — v ps [Eq. 138]. Thin solid curves lb, 2b, 3b correspond to 
impact ionization by electrons a p (E) = 0, v na — v pa [case (b), Eq. [39], dashed curves lc,2c,3c correspond to impact ionization 
by holes ce n (E) = 0, v ns = v ps [case (c), Eq. [39]. Curves of 1st, 2nd and 3rd series correspond to ao/Nd = 10~ 3 ,10~ 4 and 
1CP 5 , respectively. The parameter b/(aoEo) = 0.0002 corresponds to the doping level Nd = 10 14 cm -3 in Si. 




FIG. 6: Maximum electric field E m as a function of Vf /v s for different values of 6/ (Eocto) according to Eq. (|3Tfl . Thick solid lines 
la,2a,3a correspond to the symmetric case a n (E) — a p (E), v ns = v ps [case (a), Eq. (|38[) ]. Thin solid lines lb, 2b, 3b correspond 
to impact ionization by electrons a p (E) = 0, v ns = v pa [case (b), Eq. (|39[l ]. Curves of 1st, 2nd and 3rd series correspond to 
b/(aoEo) = 2 ■ 10 _5 ,2 • 10~ 4 and 2 • 10~ 3 , respectively; ao/Nd = 10 -4 . Insert shows the Townsends's dependence for impact 
ionization coefficient a(E). 
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FIG. 7: Slope of electric field in the screening region \d z E\ normalized by the slope in the depleted region b — qN^/eeo as a 
function of Vf/vJ for different values of ao/Nj,- Solid lines 1,2,3 correspond to the symmetric case v~ — 0. Dotted lines Id, 2d, 3d 
and dashed lines le,2e,3e correspond to the two limiting asymmetric cases v~ /vf = —1 [immobile holes, case (d)] and v~ /vf — 
+ 1 [immobile electrons, case (e)], respectively. Curves of 1st, 2nd and 3rd series correspond to ao/Nd = 10~ 3 , 10 -4 , 10~ J , 
respectively. 
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FIG. 8: Concentration of electron-hole plasma <r p i generated by the front passage as a function of front velocity Vf. In panel (a) 
the dependence (T p i(vf) is shown for different values of ao/Nd- Solid curves 1,2,3 correspond to the case of symmetric transport 
«7 = 0, fj.~ = (e.g. Vns = v ps , /J.n = (Jrp) . Dotted lines Id, 2d, 3d and dashed lines le,2e,3e correspond to the limiting cases 
of immobile holes v ps = 0, fj, p = [case(d)] and immobile electrons v n3 = 0, fj, n = [case (e)], and are calculated for the 
same values of ao/Nd- Curves of 1st, 2nd and 3rd series correspond to ao/Nd = 10 -3 , 10 -4 , 10~ 5 , respectively. In panel (b) the 
dependence cr p i(vf) is shown for different values of v^ /vf and (i~ / fi + and fixed value ao/Nd = 10~ 4 . Solid lines from 1 to 7 
correspond to v~ /vf = = —1, —0.5,0, 0.5, 0.8,0.9, 1.0, respectively. Associated dotted and dashed lines in panel (b) 

correspond to fjT //j, + = —0.9 and = 0.9, respectively, and the same value of v~ /vf as for the respective solid lines. 
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FIG. 9: Electric field _E p i in the electron-hole plasma generated by the front passage as a function of front velocity Vf. In 
panel (a) the dependence E p \(vf) is shown for different values of oo/Nd- Solid curves 1,2,3 correspond to case of symmetric 
transport vj = 0, yT = (e.g. v n3 = v ps , fj, n — fj, p ). Dotted lines Id, 2d, 3d and dashed lines le,2e,3e correspond to the 
limiting cases of immobile holes v ps — 0, n P — [case(d)] and immobile electrons v ns = 0, fj, n = [case (e)], respectively. 
Curves of 1st, 2nd and 3rd series correspond to ao/Nd = 10 -3 , 10~ 4 , 10 -5 , respectively. In panel (b) the dependence E p i(v/) 
is shown for different values of /vf and and fixed value ao/Nd = 1CP 4 . Solid lines from 1 to 7 correspond to 

v 7 / v f = A t ~// i+ — — 1) —0.5,0, 0.5, 0.8,0.9, 1.0, respectively. Associated dotted and dashed lines correspond to = —0.9 

and [f~ I fi + = 0.9, respectively, and the same value of v~ /vf as for the respective solid lines. 
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FIG. 10: Piece-wise linear approximation of the field profile used to calculate the voltage u across the n base (Sec. 3H). The 
respective cr(E) dependence is shown by dashed line A-B-C-D in Fig. [2] 



